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ABSTRACT 

Based on measurements with the Effelsberg 100-m telescope, a multi-line study of molecular species is presented toward 
the gravitational lens system PKS 1830-211, which is by far the best known target to study dense cool gas in absorption 
at intermediate redshift. Determining average radial velocities and performing Large Velocity Gradient radiative trans- 
fer calculations, the aims of this study are (1) to determine the density of the gas, (2) to constrain the temperature of 
the cosmic microwave background (CMB), and (3) to evaluate the proton-to-electron mass ratio at redshift z ~ 0.89. 
Analyzing data from six rotational HC3N transitions (this includes the J— 7^6 line, which is likely detected for the 
first time in the interstellar medium) we obtain n(H2) ~ 2600 cm~'^ for the gas density of the south-western absorp- 
tion component, assuming a background source covering factor, which is independent of frequency. With a possibly 
more realistic frequency dependence proportional to z/"'^ (the maximal exponent permitted by observational boundary 
conditions), n(Il2) ~ 1700 cm~"^. Again toward the south-western source, excitation temperatures of molecular species 
with optically thin lines and higher rotational constants are, on average, consistent with the expected temperature of 
the cosmic microwave background, Tcmb = 5.14 K. However, individually, there is a surprisingly large scatter which 
far surpasses expected uncertainties. A comparison of CS J=l<— and 4<— 3 optical depths toward the weaker north- 
western absorption component results in Tex = IIK and a 1-a error of 3K. For the main component, a comparison of 
velocities determined from ten optically thin NH3 inversion lines with those from five optically thin rotational transi- 
tions of HC3N, observed at similar frequencies, constrains potential variations of the proton-to-electron mass ratio fi 
to Afi/fj, <1.4xl0~^ with 3-(7 confidence. Also including optically thin rotational lines from other molecular species, it 
is emphasized that systematic errors are AV^ < lkms~^, corresponding to Afi/fi < 1.0x10"". 

Key words. Galaxies: abundances - Galaxies: ISM - Galaxies: individual: PKS 1830-211 - Gravitational lensing - Radio 
lines: galaxies - Elementary particles 



1. Introduction 

In the Galaxy, detailed studies of molecular clouds and 
complexes can be performed at subparsec resolution with 
a large number of spectroscopic tools in a large variety of 
regions. However, even interferometric data from nearby 
galaxies outside the Local Group only provide typical an- 
gular resolutions of a few 10 pc for quasi-thermal line emis- 
sion, matching the size of entire giant molecular clouds. 
Studies of quasi-thermal line emission from even more dis- 
tant targets are hampered by small beam filling factors 
and low flux densities. To circumvent these problems, one 
can search for masers which may be strong enough to be 
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seen out to cosmological distances (Impellizzeri et al. 2008). 
Similarly promising is the study of absorption lines. In this 
case the "effective beamwidth" is determined by the size 
of the background continuum source, which may be much 
more compact than the main beam of a single-dish tele- 
scope or the synthesized beam of an interferometer system. 
It is then the strength of the background continuum, which 
determines the sensitivity of the observations. 

In the following we will discuss spectroscopic data ob- 
tained toward a particularly prominent radio source, which 
is observed at an intermediate redshift (z~0.89), (1) to ex- 
plore the physical properties of a distant molecular cloud, 
(2) to use these properties for a test of the correlation be- 
tween Tcmb o-nd redshift^ and (3) to search for variations 
of the proton-to- electron mass ratio as a function of time. 

So far, only a few studies exist of individual giant 
molecular clouds at intermediate redshift (see, Wiklind & 



2 



Henkel et al.: Molecular absorption lines in PKS 1830-211 



Combes 1994, 1995, 1996a, 1996b; Kanekar et al. 2005). A 
good knowledge of their physical and chemical conditions 
is, however, mandatory for an evaluation of systematic dif- 
ferences between clouds at small and large cosmological dis- 
tances. 

According to the standard model of cosmology, the cos- 
mic microwave background (CMB) should follow a Planck 
curve of temperature Tcmb = (2.725±0.001) x (1-t-z), 
with < z < 108812 denoting redshift (e.g., Battistelh 
et al. 2002; Fixsen & Mather 2002; Spergel et al. 2003). 
Determining deviations from this relation would be a pow- 
erful tool to challenge the standard model of cosmol- 
ogy, highlighting possible mechanisms acting upon the 
photons of the cosmic microwave background (e.g., Lima 
et al. 2000). The best experimental method to measure 
Tcmb(1+^) is a multi-spectral line analysis of ions, atoms 
or molecules that experience negligible excitation by par- 
ticle collisions and by radiation from nearby sources, thus 
being almost exclusively excited by the isotropic radiation 
of the CMB. A comparison of line temperatures then yields 
excitation temperatures, which provide firm limits to Tcmb- 

According to the standard model of particle physics, 
fundamental constants, at least in their low-energy limits, 
should be independent of time and location. The funda- 
mental constants of physics and astronomy are well defined 
locally. Although there have been many detailed studies 
searching for spatial or temporal variations, none have been 
convincingly established (e.g., Uzan 2003; Garci'a-Berro et 
al. 2007; Flambaum 2008). However, this "constancy of con- 
stants" may not necessarily hold over the largest spatial 
and temporal scales, which are inaccessible by geological or 
astronomical studies of nearby targets. Now, advances in 
observational sensitivity make it possible to measure funda- 
mental constants for sources in the distant Universe, viewed 
at times billions of years ago and preceding the formation 
of the solar system (e.g., Bahcall et al. 1967; Varshalovich 
& Levshakov 1993; Cowie & Songaila 1995; Webb et al. 
1999). 

When searching for variations of physical constants dur- 
ing the last ~10^°yr, the most frequently studied parame- 
ter is the fine-structure constant a (Webb et al. 1999, 2001; 
Murphy et al. 2003, 2008b; Chand et al. 2004; Levshakov 
et al. 2006, 2007; Rosenband et al. 2008). Composite pa- 
rameters consisting of a, the proton-to-electron mass ra- 
tio /i, and the proton-g factor g-p have also been analyzed 
(Drinkwater et al. 1998; CariUi et al. 2000; Murphy et al. 
2001; Chengalur & Kanekar 2003; Kanekar et al. 2004, 2005; 
Tzanavaris et al. 2005, 2007). It may well be possible that 
the chromodynamic scale varies faster than the correspond- 
ing scale of quantum-electrodynamics and that, as a con- 
sequence, the proton-to-electron mass ratio ^ may show 
larger deviations from local values than the fine-structure 
constant (e.g., Calmet & Fritzsch 2002; Langacker et al. 
2002; Flambaum et al. 2004) though this is highly model 
dependent (Dent 2008). 

To date, studies focusing specifically on potential vari- 
ations of n are rare, yielding la confidence levels of order 
A/i//i ~ IQ-'^ (Levshakov et al. 2002; Ivanchik et al. 2005; 
Reinhold et al. 2006; King et al. 2008). Different powers of 
/i define the scales of electronic, vibrational, and rotational 
intervals in molecular spectra so that ro-vibronic transi- 
tions of the H2 molecule, observed toward distant quasars, 
could be used. 



An alternative and potentially even more precise ap- 
proach is the analysis of the inversion line spectrum of am- 
monia (NH3). Variations of /i could be identified as shifts 
in the relative frequencies of the inversion versus those of 
the rotational lines. After NH3 had been observed for the 
first time at significant redshift, in B0218-I-357 at z=0.685 
(Henkel et al. 2005), Flambaum & Kozlov (2007) compared 
the frequencies of the three detected inversion lines with 
those of rotational lines of CO, HCO+, and HCN (Wikhnd 
& Combes 1995; Combes & Wiklind 1995). This resulted in 
|A/i|/^ — (0.6±1.9) X 10^^. Even more recently. Murphy 
et al. (2008a) derived from a more thorough analysis of the 
ammonia spectra and newly observed HCO+ and HCN data 
A/z//x = (0.74±0.47stat±0.76sys) x 10"*^. As a cautionary 
note it has to be emphasized, however, that the NH3 inver- 
sion and HCO"*" and HCN rotational lines have frequencies, 
which differ by a factor of about ten, while the morphology 
of the background continuum is frequency dependent (e.g., 
Jethava et al. 2007). B0218-I-357 is located at a lookback 
time of 6.0 Gyr, adopting a standard A-cosmology with Hq 
= 73kms-i Mpc-\ Q.^ = 0.28, and f^A = 0.72 (Spergel et 
al. 2007). 

Among the redshifted sources showing molecular ab- 
sorption spectra, the gravitational lens PKS 1830-211 
stands out and has distinct advantages over B0218-I-357 
with respect to line strengths and chemical complexity. 
With z — 0.88582 for its main absorption feature and 
a lookback time of 7.0 Gyr, it is also the most distant 
of these sources presently known. With a lensed 2:=2.5 
quasar providing an exceptionally bright radio background, 
a vast number of molecular absorption lines could be de- 
tected (Wiklind & Combes 1996b, 1998; Gerin et al. 1997; 
Chengalur et al. 1999; Menten et al. 1999; Muller et al. 
2006). A total of 11 NH3 inversion lines has been observed 
toward the source (Henkel et al. 2008; J. Ott, priv. comm.). 
With the kinetic temperature obtained from the NH3 data, 
here we provide the complementary estimate of the density 
of the gas as well as excitation temperatures, which can 
be compared with the expected temperature of the CMB. 
To determine A/i//i, frequencies of NH3 inversion lines will 
be compared with newly measured optically thin rotational 
lines also detected in the A = 0.7-2.5 cm wavelength range. 

2. Observations 

The observations were carried out with the primary focus A 
— 1.9, 1.3, 1.0, and 0.7cm receivers of the 100-m telescope 
at Effelsberg/GermanjQ between August 2001 and March 
2002. For full width to half power (FWHP) beamwidths, 
system temperatures, and aperture efficiencies, see Table 1. 
The backend was an 8192 channel autocorrelator, which was 
split into eight segments, covering selected frequency ranges 
within a total bandwidth of 500 MHz. Each of the eight 
segments had a bandwidth of 40 MHz with 512 channels 
or, at 0.7 cm, 80 MHz with 256 channels, yielding channel 
spacings of 0.7-2.0kms~^. At 1.3 and 0.7cm, dual channel 
receivers were used. Observing with these receivers, there- 
fore at least two of the eight segments were centered at 
the same frequency to maximize sensitivity by observing 
simultaneously the two orthogonal linear polarizations. 



The 100-m telescope at Effelsberg is operated by the Max- 
Planck-Institut fiir Radioastronomie (MPIfR) on behalf of the 
Max-Planck-Gesellschaft (MPG). 
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Table 1. Observational parameters. 



A 6^ Tsys 

(GHz) (cm) (") (K) 



13.5-18.7 1.9 60-45 40 0.39-0.32 

18.0-26.0 1.3 45-35 70 0.33-0.26 

27.0- 36.7 1.0 32-23 75 0.30 

41.1- 49.7 0.7 23-19 150 0.36-0.13 



Full Width to Half Power (FWHP) beam 
width in arcsec. 

System temperatures in units of antenna 
temperature (TX). 
'^^ Aperture efficiency. 

As already mentioned (Sect. 1), PKS 1830-211 is one of 
the strongest compact continuum sources in the radio sky. 
Therefore, pointing could be checked toward the source it- 
self and was found to be accurate to about 5". The contin- 
uum was also used to calibrate the spectral lines to obtain 
line-to-continuum flux density ratios, which could be deter- 
mined with high accuracy (see Sect. 4). 

3. Results 

PKS 1830-211 is composed of a background continuum 
source at redshift 2;==2.507, possibly a blazar, and a Icnsing 
face-on spiral galaxy at z=0. 88582 along its line-of-sight 
(Wikhnd & Combes 1996b; Lidman et al. 1999; Courbin et 
al. 2002; Winn et al. 2002; de Rosa et al. 2005). The lens 
splits the background continuum into three main compo- 
nents, two compact hotspots, a north-western and a south- 
eastern one separated by '^l", and an Einstein ring, which 
is prominent at low (^10 GHz) frequencies. The bulk of 
the molecular absorption is seen toward the south-western 
hotspot (e.g., Frye et al. 1997; Swift et al. 2001). 

Figures [iHS] show the line profiles measured at Effelsberg 
with the ordinate displaying absorption in units of the ob- 
served continuum flux density. No significant indications 
for variable lineshapes were found during the time interval 
of the observations (Sect. 2), which is consistent with the 
optical depth monitoring of MuUer & Guelin (2008; their 
Fig. 5). Since about one third of the total continuum emis- 
sion is absorbed by the foreground molecular cloud at 3 mm 
wavelength (e.g., Frye et al. 1997; Wikhnd & Combes 1998; 
Swift et al. 2001), the upper panel of Fig.[l]may show a CS 
profile with intermediate (r ~ 1) peak optical depth. 

Ignoring the blueshifted feature at 14.SR~~135kms"^, 
the only spectral component known to absorb the north- 
eastern continuum hotspot (see, e.g.. Sect. 4.1.4 and 
Wikhnd & Combes 1998; MuUer et al. 2006), this CS 
J=1^0 profile may be decomposed into four velocity 
components. Gaussian fitting reveals a main feature at 
^LSR = (8±l)kms~^, a slightly weaker redshifted one 
at (18±l)kms~^, a weak narrow blueshifted feature at 
(-8±1) kms~^, and a similarly weak but broad line wing en- 
compassing several tens of kms~^. The lineshape roughly 
matches that of the opacity profiles presented by MuUer 
et al. (2006; their Fig. 5), also indicating (as the line-to- 
continuum flux density ratio) that the line is not fully 
saturated. We note, however, that the CS intensity ratio 
between the main (~8kms~^) and the redshifted compo- 



nent (~18kms~^) is slightly smaller than for the HNC and 
H^^CO+ J=:2^1 transitions; it is much smaller than for the 
HC"'^*0+ J=2^1 line, which must be optically thin. MuUer 
et al. (2006) give peak optical depths of t = 1.65, 0.27, 
and 0.15 for these three lines, respectively (their Table 4). 
The possibility that CS J=1^0 is moderately optically 
thick (t'^1) highlights the possibility that the continuum 
source covering factor may be smaller at 26 GHz than at 
3mm (see Sect. 4.1). Interestingly, the J=l<— peak line- 
to-continuum ratios of CS (Fig.[J) and C^^S (Fig.© differ 
by a factor of ~10, consistent with the corresponding ra- 
tio derived by MuUer et al. (2006; their Table 7) for the 
J=4<— 3 transition. 

The HCO+ J=1^0 profile (Fig.E lower panel), taken 
at a time when the J=2<— 1 opacity was particularly high 
(MuUer & Guelin 2008), looks quite different. Here even 
the line wings at velocities only showing weak CS absorp- 
tion are quite pronounced and the central three velocity 
components seen in CS J=1^0 appear to have merged. 
Remarkably, the velocity of the HCO+ absorption peak 
lies near -1-20 km s~^, well offset from the peak of the CS 
J=1^0 absorption. The properties of the HCO^ line may 
either be caused by a very large optical depth or by tracing 
a different volume of gas. Since the J=2^1 lines of the rare 
HCO"'" isotopologues show a "normal" optical depth profile 
(MuUer et al. 2006), a high optical depth is the preferred 
explanation. 

Figure [2] shows five detected HC3N profiles. The lines 
do not absorb more than a few percent of the total con- 
tinuum and all peak near 8kms~^. From the J=3<— 2 to 
the 6^5 transition, absolute peak line intensities rise from 
about 0.5% to 3.5% of the background continuum level, 
only to decrease further up the J-ladder. To our knowl- 
edge, it is the first time the HC3N J=7^6 line has been 
detected in the interstellar medium. With its rest frequency 
of 63.686 GHz, it is absorbed by atmospheric O2 in more 
nearby sources. 

Figure [3] shows additional spectral features which also 
absorb only a tiny fraction of the continuum background. 
These include transitions of SiO and rare isotopic species 
of CS, HCN, HCO+, and HNC. 

4. Discussion 

4.1. Excitation analysis 

Before studying molecular excitation in detail, it is impor- 
tant to obtain an estimate of the interplay between the 
radio continuum background and the intervening molecu- 
lar cloud. PKS 1830-211 is a variable source (e.g., Jin et 
al. 2003). However, for the limited time interval of our ob- 
servations we find the depths of the absorption lines and 
the continuum Uux densities to be correlated. This leaves 
the essential parameter, the line-to-continuum ratio, suffi- 
ciently unaffected to combine spectra observed at slighty 
different times (see also Wiklind & Combes 1998; MuUer et 
al. 2006; Henkel et al. 2008; MuUer & Guehn 2008). A more 
intricate problem is the background source coverage factor, 
which may vary as a function of frequency. Often, the size of 
the continuum background decreases with frequency (e.g., 
Lobanov 1998), while the size of the molecular cloud does 
not vary. If the cloud is located along the line-of-sight to 
the center of the continuum source, higher frequencies may 
yield higher source covering factors, while the opposite may 
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Fig. 1. The CS and HCO+ J=1^0 absorption lines with 
a Local Standard of Rest (LSR) velocity scale relative to 
z=0.88582, observed toward PKS 1830-211. lisR = V"hel 
+ 12.4 km s~^. Channel spacings are 0.9 and 2.0kms~^, 
respectively. The continuum level accounts for the entire 
source. Due to a highly frequency dependent noise diode 
signal at 47.294 GHz, the cahbration of the HCO+ line is 
uncertain. 



hold when the cloud is only absorbing the periphery of the 
radio background. 

At 1/ ^ 8 GHz, the absorbed south-western continuum 
source of PKS 1830-211 (see Sect. 1) contributes (38±2)% 
to the total radio flux density (e.g., Nair et al. 1993; van 
Ommen et al. 1995; Carim et al. 1998; MuUer et al. 2006; 
Menten et al. 2008). According to Frye et al. (1997), the 
molecular cloud coverage of this continuum component is 
-70%. However, Wiklind & Combes (1998) and Carilh et 
al. (1998) find values consistent with 100%, which is corrob- 
orated by the more recent monitoring results of MuUer & 
Guelin (2008). This indicates that at a frequency of 94 GHz, 
the covering factor of the total continuum is indeed fc = 
0.38±0.02. With the CS J=1^0 line, redshifted to 26 GHz 
and absorbing the continuum at a (13±1)% level (Sect. 3 
and Fig.[T]), the variation of fc between 26 and 94 GHz is 
then limited to fc oc i/O Q - 0.s_ The upper limit to the expo- 
nent would only hold, if the CS J=1^0 line were fully sat- 
urated. Because this is certainly not the case (see Sect. 3), 
the actual maximal exponent must be smaller. For t(CS 



Fig. 2. HC3N absorption lines with a Local Standard of 
Rest (LSR) velocity scale relative to z=0. 88582, observed 
toward PKS 1830-211. Channel spacings are, from top to 
bottom, 1.9, 0.7, 0.8, 1.0, 1.2, and 1.6 km s^^ The contin- 
uum level accounts for the entire source. 



J=1^0) = 1, a reasonable guess implying fc ~ 0.2 at 
26 GHz, the maximal exponent becomes 0.5. A negative 
exponent for fc as a function of frequency can be excluded 
because there is no observational evidence for larger source 
covering factors at lower frequencies (for OH and Hi, see 
Chengalur et al. 1999; Koopmans & de Bruyn 2005) and 
because then the extended, arcsecond sized Einstein ring 
becomes prominent (e.g., Jauncey et al. 1991). 

With five detected and one undetected rotational tran- 
sition, we can carry out the first multilevel HC3N study of a 
significantly redshifted source (see Mauersberger et al. 1990 
for an analogous study of a nearby starburst galaxy). With 
its small rotational constant, HC3N is particularly useful 
to constrain molecular excitation along the line-of-sight to 
PKS 1830-211. Two SiO transitions and comparisons of our 
data with spectra observed at higher frequencies provide 
additional constraints. The analysis is based on measured 
"apparent" optical depths. 



-In 1 



Tc 



with Tl and Tq denoting line and continuum temperature 



and with t. 



app 



\T-l\/Tc being a good approximation in 
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Table 3. HC3N peak optical depths corrected for incom- 
plete continuum source coverage (see Sect. 4.1.1). is the 
source covering factor, v denotes the redshifted frequency. 
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Fig. 3. Additional optically thin absorption lines with a 
Local Standard of Rest (LSR) velocity scale relative to 
z=0.88582, observed toward PKS 1830-211. Channel spac- 
ings are, from top to bottom, 1.0, 2.0, 0.9, 2.0, 2.0, and 
2.0kms^^. The adopted continuum level includes the en- 
tire source. 



the case of optically thin lines (beware of the typo in the 
corresponding equation of Henkel et al. 2008). 



Table 2. Line widths and optical depths of HC3N transi- 
tions (see Sect. 4.1. 1)''^ 



Line 






FWHP 


7~pcak 


\'tAV 


HC3N J= 


=3* 


-2 


10.6±1.2 


0.0057±0.0008 


0.060±0.005 


HC3N J= 


=4* 


-3 


9.3±1.0 


0.0131±0.0018 


0.121±0.010 


HC3N J= 


=5^ 


-4 


8.6±1.0 


0.0211±0.0031 


0.181±0.017 


HCaN J= 


=6^ 


-5 


9.0±1.9 


0.0280±0.0069 


0.251±0.034 


HC3N J= 


=7^6 


12.2±2.1 


0.0173±0.0035 


0.210±0.024 


HC3N J= 


10 


^9 






<0.235 



Line 


/c = 0.3 


/c= 0.2x(!//26GHz)"=' 


HC3N J=3^2 


0.0239±0.0046 


0.0462±0.0089 


HC3N J=4^3 


0.0482±0.0093 


0.0807±0.0156 


HC3N J=5^4 


0.0721±0.0106 


0.1084±0.0158 


HC3N J=6^5 


0.1000±0.0221 


0.1367±0.0302 


HC3N J=7^6 


0.0837±0.0174 


0.1059±0.0219 



''^ Full Width to Half Power (FWHP) line widths, peak appar- 
ent optical depths, and integrated optical depths obtained 
from one component Gaussian fits. The upper limit to the 
HC3N J=10<— 9 integrated opacity is a 3o" value. 



4.1.1. HC3N 

To reproduce the HC3N data (Fig. [2]), we use a Large 
Velocity Gradient (LVG) model including 22 levels up to 
92 K above the ground state. Collision rates were taken 
from Green & Chapman (1978). The bulk of the ammo- 
nia (NH3) column density arises from gas at Tkin '--^ 80K 
(Henkel et al. 2008). This value was adopted. As already 
mentioned (Sect. 3), the HC3N lines only absorb a few per- 
cent of the total continuum background (Fig. [51 Table 2). 
This is similar to the NH3 lines also seen at cm- wavelengths. 
Since the NH3 lines are definitely optically thin (the hy- 
perfine satellites are at best barely visible), since HC3N 
and NH3 frequencies are similar, and since the continuum 
source covering factor is of order /c = 0.10-0.38 (Fig. [T] and 
Sect. 4.1), we conclude that all detected HC3N lines must 
also be optically thin. 

For the LVG code, a spherical cloud geometry was as- 
sumed, but with the lines being optically thin, results do 
not depend on cloud morphology. We consider the inte- 
grated optical depths to be more reliable than the peak 
optical depths. For the input parameters we therefore only 
adopted the peak optical depth of the very well fitted 
J=5<— 4 line, accounting for a continuum source covering 
factor /c = 0.3 and using the ratios of the integrated optical 
depths to obtain the peak optical depths of the other lines. 
This results, ignoring potential uncertainties in the contin- 
uum coverage, which are difficult to quantify, in the peak 
optical depths displayed in the central column of Table 3. 

As a first approach we neglect coUisional excitation 
and only vary the temperature of the cosmic microwave 
background to simulate the observational data with min- 
imized reduced values, Xred = X^/{^~P) {N=5: num- 
ber of lines, P—2: number of free parameters; see Jethava 
et al. 2007). At a redshift of z=0.88582, we expect Tcmb 
= (2.725±0.001) X {l+z) = (5.139±0.002)K. With this 
value, the minimal x^^^ is an unacceptable 8.1 (Table 4). 
Modifying Tcmb, but still only permitting excitation by 
the cosmic microwave background, the optimal value is ob- 
tained with Xred = 0-8 at Tcmb ~ 20K (model A7). Such 
a high Tcmb value is also unacceptable. We therefore have 
to conclude that our model assumptions are too simple. 

As the next step to higher complexity, we therefore also 
include collisions with H2. Keeping Tcmb fixed at 5.14, we 
calculate the radiative transfer for 500 cm~^ < n(H2) < 
5000 cm~^. Already at n(H2) = 500 cm"^, the quality of the 
fit is clearly improved (model B2 of Table 4). The optimum 
is reached at n(H2) ~ 2600 cm^'^ (models B7/B8) with x^cd 
= 0.4. In this case, calculated and measured optical depths 
differ by about 15%. Since the opacity ratios do not vary 
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Table 4. 


HC3N 


LVG models 


and 


resulting 


reduced x 


^ values (see Sect. 4.1.1)°. 




Model 


-tCMB 


n(H2) 


Xrcd 












(K) 


(cm-^*) 






Al 


5.14 





8.1 


4.1 


A2 


6.00 





6.0 


2.4 


A3 


8.00 





3.3 


0.8 


A4 


10.00 





2.0 


0.5 




12.00 





1.4 


0.6 


A6 


15.00 





1.0 


1.0 


A7 


20.00 





0.8 


1.7 


A8 


25.00 





0.9 


2.2 


A9 


30.00 





1.0 


2.6 


AlO 


50.00 





1.4 


3.6 


Bl 


5.14 





8.1 


4.1 


B2 


0.14 


500 


5.4 


2.1 


B3 


5.14 


1000 


3.3 


0.8 


B4 


5.14 


1500 


1.7 


0.3 


B5 


5.14 


1700 


1.0 


0.2 


B6 


5.14 


2000 


0.8 


0.3 


B7 


5.14 


2500 


0.4 


0.7 


B8 


5.14 


2700 


0.4 


1.0 


B9 


5.14 


3000 


0.5 


1.6 


BIO 


5.14 


4000 


2.1 


4.0 


Bll 


5.14 


5000 


4.9 


7.0 



'Assumed kinetic temperature: Tkin = 80 K. 
For the reduced values with the frequency 
independent continuum source coverage factor 
fc ~ 0.3, see Col. 4. Accounting for the max- 
imal frequency dependence of fc, fc = 0.2 x 
(!^/26 GHz)"'^, the corresponding reduced 
values are given in the last column. A com- 
parison of the results summarized in the last 
two columns provides a good measure of the 
uncertainties involved. 



with column density as long as the lines are optically thin, 
this result is not sensitive to the specific choice of fc {fc oc 
r-i cx N^Csn)- Already at n(H2) = SOOcm^^ the HC3N 
J=1^0 line is inverted. At ^(Hs) = 2000 cm'^, the 2^1 
line is also inverted and at n(H2) = 7000 cm~'^, the three 
ground rotational lines of HC3N are all inverted. Thus, at 
n{ll2) ^ 7000 cm"'^, we do not expect to obtain a detectable 
absorption signal in the J— 3^2 transition. 

Instead of applying a constant fc — 0.3, we can also use 
the opposite extreme, an fc value with the strongest fre- 
quency dependence permitted by the boundary conditions 
discussed in Sect. 4.1. For fc = 0.2 x (i//26 GHz)^-^ the 
resulting optical depths are given in the right column of 
Table 3. With fc decreasing with decreasing frequency, the 
lower HC3N transitions become more prominent (Tj-cai oc 
/~^) so that less excitation is required to simulate the calcu- 
lated optical depths. The right column of Table 4 indicates 
this change. Excluding collisions, the lowest reduced is 
reached at Tcmb = lOK instead of the previously obtained 
20 K. With Tcmb = 5.14K and introducing collisions, we 
get n(H2) ~ 1700 cm~'^ instead of 2600 cm^'^. Again, cal- 
culated and measured hue opacities differ by about ±15%. 
Overall, the Xrcd values tend to be smaller than those ob- 



tained with fc = 0.3. This may indicate that the approach 
with variable fc is the better one. 

To compare our results with those for the slightly less 
redshifted source B0218-h357 (Henkel et al. 2005, 2008; 
Jethava et al. 2007), we conclude that in PKS 1830-211 
kinetic temperature and gas density are much higher, while 
the source covering factor may be less frequency dependent. 
However, the situation may not be quite as clear as it ap- 
pears at first sight. The kinetic temperatures were mea- 
sured with the same tool, namely the inversion lines of am- 
monia, thus providing an objective comparison. The gas 
densities, however, are determined in different ways, ana- 
lyzing formaldehyde (H2CO) in the case of B0218-I-357 and 
cyanoacetylene (HC3N) in the case of PKS 1830-211. This 
may lead to a systematic bias. If we were able to detect 
HC3N absorption in its ground rotational transition toward 
PKS 1830-211, for example, the density of the gas traced by 
this line should only be a few 100 cm~'^ and not 2000 cm^"^. 
Otherwise the level populations would be inverted. While 
HC3N appears to be very weak in B0218-I-357, an H2CO 
study hke that of Jethava et al. (2007) for B0218-I-357 is 
also feasible for PKS 1830-211 to provide a comparison of 
densities excluding any potential bias caused by the use of 
different molecular tracers. 

In view of our goal to constrain Tcmb (see Sects. 4.1.2 
and 4.1.3), it is important that the density derived is not 
sensitively dependent on the choice of the temperature of 
the cosmic microwave background. To test this, we also cal- 
culated models with Tcmb = 4K and 7K instead of 5.14 K. 
Resulting densities are n(H2) ^ 2000 cm~^ and ~1000cm~'^ 
instead of ~1700cm~'^, when using frequency dependent 
source covering factors. 

4.1.2. SiO 

HC3N was particularly interesting, because it is a heavy 
molecule with a small rotational contant. In PKS 1830-211, 
HC3N rotational lines can be observed at multiples of ~ 
4.8 GHz, allowing us to detect several transitions. Small fre- 
quencies lead to slow rates of spontaneous decay (Ay oc v^) 
and level populations are already altered by coUisional pro- 
cesses in a rather tenuous interstellar medium with n(H2) 
- 1000cm"3. 

There exist a large number of lighter interstellar 
molecules with higher rotational constants and low appar- 
ent optical depths, indicating as in the case of HC3N op- 
tically thin absorption. SiO is one of these species, which 
exemplifies some effects we have to account for in a thor- 
ough analysis of PKS 1830-211. SiO has a five times larger 
rotational constant than HC3N and a similar electric dipole 
moment. With Tcmb — 5.14 K and a density of 2000 cm~'^ 
(Sect. 4.1 and Table 4), LVG modeling (collision rates from 
Turner et al. 1992) yield in the optically thin limit excita- 
tion temperatures of 5.51 K and 5.21 K for the SiO J=1^0 
and 2<— 1 lines, respectively. This implies that the densities 
derived in Sect. 4.1.1 are not drastically modifying line in- 
tensities and that, if SiO arises from the same gas as HC3N, 
the bulk of the excitation is provided by the CMB. 

Line-to-continuum ratios of 0.13 or 0.38 are observed 
in the CS 1^0 (Fig.[l|) or HON and HCO+ J=2^1 lines 
(MuUer & Guelin 2008), providing firm boundary condi- 
tions for the frequency range from the SiO 1<— to the 
2<— 1 transition. Since the SiO lines are much weaker, they 
must be optically thin. With integrated apparent opti- 
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Table 5. Integrated optical depths of optically thin lines (see Sects. 4.1.2 and 4.1.3) 
and resulting excitation temperatures. 







JV dV 




Excitation temperature 




J=1^0 


7=2^1" 


J=2«-l/J=1^0'" 


rph 


-'cx,2-l 


^ cx.LTE 




(kms 








(K) 




SiO 


0.412±0.021 


1.440±0.270 


3.495±0.679 


5.51 


5.21 


23 




0.292±0.020 


0.440±0.040 


1.507±0.172 


5.83 


5.27 


7.2±0.4 


H"CO + 


2.804±0.259 


2.500±0.050 


0.892±0.084 


5.20 


5.17 


3.8±0.3 


H"CN 


1.791±0.243 


2.100±0.080 


1.173±0.165 


5.37 


5.24 


4.8±0.5 


HN"C 


0.307±0.181 


0.900±0.070 


2.932±1.743 


5.20 


5.17 




CS-NE'' 


0.080±0.010 


0.280±0.080 


3.500±1.092 






11±3 



^-^ In the cases of C'^'^S and CS, J=2^1 should be replaced by J=4^3. With the exception 
of the SiO J=2<— 1 line (Fig.[3ll, the J=2^1 and 4<— 3 data were taken from MuUer et al. 
(2006). 



LVG excitation temperatures of optically thin lines with Tkin = 80 K, Tcmb = 5.14 K, and 
n(H2) = 2000 cm-^ 

Excitation temperatures under conditions of Local Thermodynamical Equilibrium (LTE). 
For the uncertainty in the estimate from SiO (this includes so far unpublished SiO J=4-3 
data), see Sect. 4.1.2. 

' CS absorption toward the north-eastern continuum component, displaced from the main 
absorption line by -147kms~^. Since no spatial density has been derived for this feature. 
Cols. 5 and 6 are left empty. 



cal depths of (0.412±0.021)kms-i and (1.44±0.27) kms^^ 
for the 1^0 and 2^1 lines, respectively, we obtain 
r(2^1)/r(1^0) = 3.50±0.68 (Table 5). Assuming Local 
Thermodynamical Equilibrium (LTE), this results in an 
excitation temperature of 23 K, with the expected 5 K (re- 
quiring a ratio of 2.2) at a 2a level. We conclude that the 
noise in the SiO J=2^1 line is still too high for a definite 
result. More sensitive J=2^1 data are therefore required. 

Using SiO as an example, we should also note that there 
exists the possibility that the absorption exclusively arises 
from a much denser gas component than HC3N. Unlike 
emission lines, the strength of absorption lines does not 
depend on "critical densities" , where timescales of radia- 
tive de-excitation match those for coUisional excitation. 
Therefore, a low excitation level caused by large frequen- 
cies and dipole moments does not necessarily lead to non- 
detections and a high ratio of the optical depths between 
the 2*-l and l<-0 lines of SiO at Tcmb ~ 5.14 K could be 
the consequence of chemical differentiation. A test with a 
more recently taken SiO J=4-3 line [j rdv = 0.70±0.07; 
S. MuUer, priv. comm.) yields Trot = (6.8±0.3)K. While 
this may indicate that there is no dense SiO component, 
we emphasize that half a decade lies between this recent 
measurement from August /September 2006 and the other 
data presented here. Variability of the source might thus 
lead to a larger uncertainty than the given formal error. 

Tests with our LVG model show another effect. The ex- 
citation of the J=l<— line rises with increasing density, 
becoming inverted at n(H2) ~ 6xl0^cm~'^. While rising, 
the J=l<— optical depth decreases, keeping the product 
rTex approximately constant. Thus r(2<— l)/r(l<— 0) ratios 
in excess of four (the maximum for LTE conditions) are 
possible and may become a sensitive probe of densities re- 
lated to molecular constituents, unaffected by critical den- 
sities. For Tkin = 80K and Tcmb = 5.14K, a density of 
1.3x 10** cm^'^ would reproduce the measured but uncertain 



ratio of '^3.5. This density is almost an order of magnitude 
higher than that derived from HC3N. 

4.1.3. Other molecular species 

Here we will analyze opacity ratios between optically thin 
J=l<— lines measured by us (Fig. [3]) and the corresponding 
J=2^1 or J=4^3 lines previously reported by MuUer et 
al. (2006; their Table 4). CS has a rotational constant sim- 
ilar to that of SiO. HCN, HCO+, and HNC have rotational 
constants about twice as large and an order of magnitude 
larger than that of HC3N. This drastically reduces effects 
caused by coUisional excitation as long as these species arise 
from the same volume as HC3N. Columns 5 and 6 of Table 5 
show this effect. With the boundary conditions of Tkin = 
80 K, Tcmb = 5.14 K, and ^(Hs) = 2000 cm'^ (Sect. 4.1.1), 
excitation temperatures for HCN, HCO+, and HNC are 
only slightly larger than Tcmb, while for CS the difference 
is a little larger. The values were obtained from LVG calcu- 
lations with collision rates of Turner et al. (1992), Flower 
(1999), and Schoier et al. (2005). 

While the signal-to-noise ratio of our HN^^C J=l<— 
line is not high enough to derive a meaningful value, 
the average of the LTE excitation temperatures of C'^'^S, 
H^^CO+, and H^^CN is close to the expected value of 
5.14 K. However, deviations of individual values (see Col. 7 
of Table 5) are significant in most cases and are well beyond 
the estimated uncertainties. Combes & Wiklind (1999) 
showed a diagram of previously determined LTE temper- 
atures (their Fig. 4), obtained with the SEST 15-m and 
IRAM 30-m telescopes. These results also show significant 
scatter. While LTE temperatures higher then the expected 
Tcmb value could be explained by additional coUisional or 
radiative excitation (in the case of our C'^^'S results, a den- 
sity of n(H2) = 10"* cm^'^ would be needed), those ratios 
below the expected Tcmb value, obtained by Combes & 
Wiklind (1999) and by us from linear molecules, are more 
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difficult to interpret. Combes & Wiklind (1999) suggested 
that their low values, not including J=l<— transitions, are 
caused by the heating of the 1^0 and an associated cool- 
ing of the 1 transition. Our data do not conffi'm this 
hypothesis, indicating instead that the scatter in the LTE 
excitation temperatures is much higher than the errors in 
Col. 7 of Table 5 suggest. A dedicated effort including many 
more species, among them also asymmetric tops, may solve 
this problem. So far, the molecular absorption lines toward 
PKS 1830-211 do not provide the kind of constraints on 
TcMB, which one is also attempting to obtain in analyzing 
atomic fine-structure levels in redshifted sources (e.g., Ge 
et al. 1997; Srianand et al. 2000; Molaro et al. 2002). 

4.1.4. The northwestern source 

As already pointed out in Sect. 3, absorption toward 
the south-western compact continuum component of 
PKS 1830-211 is relatively strong. However, a few promi- 
nent molecular lines are also detected at a velocity offset 
of -147kms~^ toward the north-eastern component, which 
is located behind another spiral arm of the lensing galaxy 
(e.g., Wiklind & Combes 1998; MuUer et al. 2006). While 
our HCO"'" profile (Fig.[T|) is not calibrated well enough for 
an analysis, the CS J=1^0 feature (also Fig.[T]) can be 
compared with the corresponding J=4<— 3 profile of MuUer 
et al. (2006). The resulting LTE excitation temperature is 
given in the last line of Table 5. As a caveat, however, we 
note that the CS spectrum from MuUer et al. (2006) was 
taken in July 1999, more then two years before our mea- 
surements (see MuUer & Guelin 2008 for source variability 
with time) . This time interval is longer than for the spectra 
from the south-western source. 



4.2. The proton-to-electron mass ratio 
4.2.1. Some general considerations 

As already mentioned in Sect. 1, a comparison of radial ve- 
locities from inversion lines of ammonia (NH3) with rota- 
tional lines of either NH3 or other species can provide sensi- 
tive constraints on variations in the proton-to-electron mass 
ratio fj, = m-p/mc- Following Flambaum & Kozlov (2007), 



c 



l + z 



(1) 



with zjnv and Zjot denoting the apparent redshifts of the 
inversion and rotational lines, z representing the real red- 
shift of the source, and A/z providing the deviation relative 
to the current laboratory value of /z, defined such that a 
negative value indicates a smaller in the absorbing cloud. 
In the case where both the NH3 and rotational transitions 
show absorption in just a single velocity component, defin- 
ing and measuring the velocity difference IS.V is straight- 
forward. The more general case where many velocity com- 
ponents exist must be treated differently, by varying all the 
free parameters, including AT^, to minimize in a simul- 
taneous fit to all NH3 and rotational profiles (e.g.. Murphy 
et al. 2008a). 

Ammonia and cyanoacetylene (HC3N) are both molec- 
ular high density tracers. Toward PKS 1830-211, the NH3 
inversion lines and the rotational HC3N lines are optically 
thin and observed at similar frequencies. Therefore, in this 
paper we derive A/i//i by comparing the NH3 and HC3N 



Table 6. Radial velocities (see Sect. 4.2). 
Uncertainties in the frequencies are small relative 
to those of the velocities, which were obtained by 
Gaussian fits". 



Species 



Transition 



Rest 
Frequency 
(GHz) 



Vlsr 
(kms-i) 



NH3 



{J,K) 



(1,1) 
(2,2) 

(3,3) 
(4,4) 

(5,5) 

(6,6) 

(7,7) 
(8,8) 
(9,9) 
(10,10) 



HC3N 



7^6 



SiO 
SiO 



23.694496 

23.722633 

23.870130 
24.139416 

24.532989 

25.056025 

25.715182 
26.518981 
27.477943 
28.604737 

27.294289 
36.392324 
45.490314 
54.588247 
63.686052 

43.423846 
86.847010 
48.206946 
86.754288 
86.339922 
87.090873 



9.0 
8.6 
7.4 

10.5 
8.6 
8.4 
9.1 
8.2 
8.7 
8.1 
8.5 

10.9 
8.5 
8.8 

10.2 

9.0 
8.2 
8.0 
7.7 
10.0 

9.8 
8.9 
9.8 
9.7 
8.8 
7.4 



0.1 
0.5 
0.2 
0.4 
0.1 
0.2 
0.7 
0.3 
0.5 
0.2 
0.6 
0.4 
0.6 
0.7 
2.0 

0.4 
0.3 
1.0 
0.3 
0.4; 

0.7 
1 

0.6 
0.7 
1.5 
1.3 



Frequencies were taken from the JPL cata- 
log (NH3), the Cologne Database for Molecular 
Spectroscopy (HC3N), and the line list of Lovas 
(other species). NH3 transitions with two Vlsr val- 
ues were measured twice (see Henkel et al. 2008). 



spectra. Given the fairly low signal-to-noise ratio of these 
spectra, there is little obvious evidence for more than a 
single absorbing velocity component across the NH3 and 
HC3N profiles. We therefore proceed with a simple analy- 
sis assuming that only one absorbing velocity component 
exists, deferring a simultaneous minimization analysis 
which incorporates other weak components and other rota- 
tional transitions to a later paper. Despite this assumption 
of a single velocity component, judicious treatment of the 
uncertainties on its velocity in the NH3 and rotational tran- 
sitions will ensure our constraint on A/i//i is robust, if not 
optimally precise. 

The NH3 inversion lines were measured with the Green 
Bank Telescope (GET), while the HC3N data were obtained 
at Effelsberg. The Local Standard of Rest velocity scales at 
both sites are, however, well adjusted (Tifft & Huchtmeier 
1990) and have not been significantly modified in recent 
years. They refer to a peculiar solar velocity of 20kms~-'^ 
toward aigoo = 18*^ and 

(51900 = +30°, corresponding to 
14sR - ^HEL = 12.432 km s"i in the case of PKS 1830-211. 
Therefore, the velocity measurements listed in Table 6 can 
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be directly compared to estimate A/i//i and its uncertainty, 
without further velocity corrections. 

4.2.2. Average velocities of NH3 and HC3N 

For the single-Gaussian fits to the NH3 transitions, Table 6 
reveals somewhat inconsistent centroid velocities. They 
scatter around the weighted mean velocity over all tran- 
sitions of 8.65 km s^^ with a value of 117 which, for 
just 14 degrees of freedom, is far more than expected based 
on the individual velocity uncertainties. This may indi- 
cate the presence of more than one velocity component. 
Despite the large range of NH3 excitation states observed 
by Henkel et al. (2008) - some 1030 K - we do not observe 
a trend in velocity with excitation state: For the three low- 
est NH3 inversion lines, the unweighted mean velocity is 
8.82±0.45kms~^ and for the higher excited ones we find 
a mean of 8.94±0.27kms~^ (the errors are standard de- 
viations on the means). Indeed, the larger than expected 
scatter in fitted velocities for the different NH3 transitions 
seems to be random in sign and magnitude. A similar ad- 
ditional scatter is also observed in the rotational HC3N 
transitions (Fig. [2]). These cover only a small amount of 
excitation, with the highest level, J=7, being just ~12K 
above the ground state. 

Given the velocity scatter in the transitions of NH3 and 
HC3N, constraining the variation of /i is best done with 
conservative velocity and uncertainty estimates which nat- 
urally incorporate the observed scatter. The simplest such 
velocity estimator is obviously the unweighted mean veloc- 
ity and its standard deviation which, for the NH3 transi- 
tions in Table 6, are 8.90±0.24kms"\ and for the HC3N 
transitions, are 8.58±0.37kms~^. Of course, the statistical 
velocity uncertainties quoted in Table 6 for some transi- 
tions are so high that including them in the mean velocity 
calculation is likely to decrease the reliability of the mean. 
We therefore reject transitions with velocity uncertainties 
larger than the root-mean-square (RMS) velocity variation 
for each species, 0.91 and 0.83 km s^^ for NH3 and HC3N re- 
spectively. With this criterion, the NH3 (10,10) and HC3N 
5<— 4 transitions are rejected. The final mean velocity for the 
remaining NH3 transitions is 8.81zL0.23 km s~^ and, for the 
remaining HC3N transitions, is 8.73zL0.43 kms~^. 

4.2.3. Resulting limits for a variation of ^ 

Using these final clipped mean velocities, AV — 
0.08±0.49kms^^. Equation [1] then provides our l-cr con- 
straint on the variation in fj,, A^/^ = (-|-0.08±0.47) x 10^^. 
Since the quoted uncertainty derives entirely from the scat- 
ter in the individual transition velocities, this should be a 
reasonably robust error estimate. Nevertheless, given that 
we have used only single Gaussian fits to the absorption 
profiles and that there is some scatter in the individual 
transition velocities, we quote our final result as a 3-cr up- 
per limit on variation in /i at the absorption redshift of 
z = 0.89, 

lA^i/^il < 1.4 X 10"^ (2) 

For comparison with laboratory constraints on variations in 
yU, and in the absence of a reliable model for how fi might 
be expected to vary with cosmological time, it is common, 
if not well motivated, to assume that any variation is linear 



in time. Hence, our upper limit on variation in fi translates 
to a 3-cr upper limit on its time variation 

< 2.0 X IQ-^^yi-^ (3) 

over the past 7.0 Gyr. 

There is yet another study on the proton-to-electron 
mass ratio in the main lens of PKS 1830-211. From a com- 
parison of the ammonia inversion lines (Henkel et al. 2008) 
with its { J,K) ~ (1,0)<— (0,0) rotational transition, Menten 
et al. (2008) find consistency, at a la level, of Afi/fi = 
1.9xl0~^. The strength of this study is its focus on lines 
with different dependencies on fj, arising from the same 
molecular species. On the other hand, the ratio between the 
frequencies of the rotational and inversion lines is '^25, thus 
leading to potentially significant differences in the mor- 
phology and the covering factor of the background radio 
and submillimeter continuum. While such systematic dif- 
ferences cannot be quantified on the basis of a single ro- 
tational line, the given uncertainty of the resulting A/Lt//Lt 
value is dominated by the limited signal-to-noise ratio of 
the rotational line. 

4.2.4. Potential caveats 

The consistency of NH3 and HC3N radial velocities is re- 
markable in view of a number of effects, which might ex- 
ert a significant influence on our results. Having carefully 
avoided the use of optically thick transitions, these are (1) 
time variability of the continuum source, (2) a frequency 
dependent continuum morphology, (3) hyperfine structure, 
(4) chemistry, and (5) inhomogeneities in temperature and 
density inside a cloud of size ^10 pc (Carilli et al. 1998). 

(1) A time variable continuum source may lead to dif- 
ferent lines-of-sight and thus to different radial velocities 
at different epochs. Not much of this was seen in spite of 
the existence of monitoring programs (Wiklind & Combes 
1998; MuUer et al. 2006; MuUer & Guehn 2008) or repeated 
measurements of ammonia lines (Table 6 and Henkel et al. 
2008). 

(2) The NH3 and HC3N transitions used in this study 
are much closer in frequency than those chosen by Murphy 
et al. (2008a) and Menten et al. (2008). Also, in PKS 1830- 
211, the source covering factor may be less frequency depen- 
dent than in B0218-t-357 (Sect. 4.1.1). Nevertheless, NH3 
and HC3N frequencies differ by factors of 1.0-2.5. While 
no direct effect is apparent, it remains a source of uncer- 
tainty which cannot be quantified. 

(3) In particular NH3 velocities would be greatly af- 
fected by hyperfine (hf) splitting, if strong non-LTE ef- 
fects would occur. A main group of hf-components is sur- 
rounded by four satellite groups, displaced by about ±10 
and ±20kms^^. The higher the energy above the ground 
state of an inversion line, the weaker the satellite features 
relative to the main group (e.g., the satellites account for 
^2.5% of the total absorption in the {J,K) = (7,7) line). 
The above made comparison of average velocities of the 
three lowest with the seven higher excited NH3 inversion 
lines does not show a significant shift in velocity. Deviations 
from LTE intensity ratios of the various components are 
only expected in the case of a significant population of the 
non-metastable inversion states {J>K), which requires ex- 
tremely high densities or intense radiation fields (Stutzki 
et al. 1984; Stutzki & Winnewisser 1985a,b). Non-LTE ef- 
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fects leading to HC3N velocity shifts as a function of ro- 
tational quantum number J are also not obvious. In this 
case, the J=3^2 line would be the most critical, while the 
hf-structure of the T-i— 6 line is far too compact to yield any 
significant shifts. 

(4) Chemistry is a more difficult issue, because there 
is no way to estimate its effect on our data. NH3 and 
HC3N may have different spatial distributions, which is 
most impressively exemplified by the nearby prototypical 
dark cloud TMC-1, where peaks of line emission are offset 
by Tarcmin (e.g., Olano et al. 1988). HC3N is a molecule 
representing a young chemical age, while NH3 refers to 
a chemically more evolved stage of cloud chemistry (e.g., 
Dickens et al. 2001; Hirota et al. 2002). 

(5) Since NH3 is covering an enormous range of excita- 
tion without showing any significant change in radial veloc- 
ity, changes in temperature and density may also not play 
an important role. 

4.2.5. Summary and outlook 

To summarize, we either have the unlikely scenario in which 
different effects cancel each other out, or all of these effects 
are minor. In order to broaden the scope of our study and 
to avoid relying entirely on the rotational transitions of 
HC3N, we may also use the J=l<— and 2^1 lines pre- 
sented in Fig. [31 We should note, however, that for these 
lines, frequencies are significantly higher than for NII3 and 
HC3N. The unweighted mean rotational velocity then be- 
comes 9.07zL0.38 km . The velocity difference amounts to 
AT/ = +0.26±0.U kms~^^ relative to NH3, while for HC3N 
we have found -0.08±0.49 km . These differences high- 
light potential systematic effects, caused by a time variable 
continuum background, by a frequency dependent contin- 
uum morphology, by non-LTE effects in lines with hyper- 
fine structure, by astrochemistry, and by the specific exci- 
tation requirements of individual lines. The large number of 
analyzed molecular transitions, covering a wide parameter 
space, ensures that systematic discrepancies in velocity are 
within Ikms""'^, corresponding to A///// <10~^. 

Our constraint on variations of the proton-to-electron 
mass ratio /i differs in several ways from that given by 
Murphy et al. (2008a) for B0218-I-357 at z~0.68 (this study 
was also based on a comparison of ammonia (NH3) in- 
version lines with rotational transitions of other molecular 
species). (1) The frequencies of the lines measured by us are 
much more similar than in the previous study, minimizing 
effects caused by a frequency dependent continuum back- 
ground morphology on the observed column. (2) All our 
analyzed lines are optically thin, avoiding modified line- 
shapes due to saturation and (3) our number of observed 
transitions is much larger. This may result in an evalu- 
ation of hidden systematic errors, which is more realistic 
than any estimate from multi-velocity component fits of a 
few spectral lines only. (4) The time interval between our 
measurements of inversion and rotational lines is smaller, 
possibly reducing systematic effects caused by a time vari- 
able radio continuum background. Note, however, that this 
may be compensated by a higher variability of our source. 

There are also two weaknesses with respect to the study 
of Murphy et al. (2008a). (1) The excitation of most inver- 
sion lines analyzed by us is much higher than that of the 
rotational lines. While no trend of radial velocities as a 
function of excitation is seen by us, the study of Murphy et 



al. (2008a) is confined to inversion lines of modest excita- 
tion and is definitely less sensitive to potential discrepan- 
cies with respect to excitation. (2) We did not yet provide a 
multi-component spectral analysis of the entire dataset, ac- 
counting for velocity structure and hyperfine splitting. This 
is the missing step toward a complete analysis of the spec- 
tra and will be discussed in a forthcoming paper, focusing 
entirely on potential variations of 11. 

5. Conclusions 

Our observations of molecular lines at redshift Z'^0.89 to- 
ward the gravitational lens system PKS 1830-211 reveal the 
following main results: 

— Cyanoacetylene (IIC3N) is detected in absorption in five 
rotational lines. The J=7<— 6 line may have been de- 
tected for the first time in the interstellar medium. CS, 
C^^S, HCO+, HC"0+, Hi^CN, and SiO absorption are 
also detected, in their ground rotational transitions. 

— A Large Velocity Gradient analysis of the HC3N spec- 
tra yields, adopting Tkin = 80 K from ammonia (NH3), a 
density of n(H2) ^ 2600 cm^"^ for a background source 
covering factor /c, which is independent of frequency. 
CS, HCO"*", and OH data constrain its frequency de- 
pendence to /c oc i/0 - 0.5^ This holds for the frequency 
range between 26 GHz and 94 GHz. For a maximal fc 
exponent of 0.5, n(H2) is reduced to ~1700cm~'^. 

— The gas density determined using HC3N is low enough 
that molecules with higher rotational constants should 
be mainly excited by the cosmic microwave background. 
Therefore, we compare our C^^S, R^^CO+, and H^^CN 
J=1^0 data with the corresponding J=2<— 1 lines 
( J=4^3 in the case of C'^^'S) presented by MuUer et al. 
(2006). While the average excitation temperature is con- 
sistent with the expected temperature of the microwave 
backgrund, Tcmb = 5.14K, deviations of individual val- 
ues (3-7 K) far exceed expected uncertainties. 

— For the weak north-eastern absorption component, a 
comparison of the CS J=l<— and 4<— 3 transitions 
yields an excitation temperature of (11±3) K. The com- 
parison is, however, based on data, which were taken at 
epochs more than two years apart. 

— From a comparison of the radial velocities of five rota- 
tional transitions of HC3N with ten inversion transitions 
of ammonia, the proton-to-electron mass ratio n of the 
lens at z=0.89 deviates by |A^/^| < 1.4 x 10~^ from 
the local value with 3-cr (99.7%) confidence. Assuming 
that any variation in fj, evolves linearly with cosmolog- 
ical time, this corresponds to < 2.0 x 10^^^ yr^^. 
The large number of analyzed optically thin molecu- 
lar transitions, covering a wide range of physical and 
chemical parameters, indicates that systematic effects 
should be <lkms~^, corresponding to an uncertainty 
of < 10~s in A/i//i. 

— Comparing the main lens of PKS 1830-211 with that of 
B0218+357, the molecular gas toward PKS 1830-211 is 
warmer and appears to be denser, while its molecular 
source covering factor may be slightly less frequency 
dependent at cm- and mm-wavelengths. 
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